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We use a numerical implementation of the strong disorder renormalization group (RG) method 
to study the low-energy fixed points of random Heisenberg and tight-binding models on different 
types of fractal lattices. For the Heisenberg model new types of infinite disorder and strong disorder 
fixed points are found. For the tight-binding model we add an orbital magnetic field and use both 
diagonal and off-diagonal disorder. For this model besides the gap spectra we study also the fraction 
of frozen sites, the correlation function, the persistent current and the two-terminal current. The 
lattices with an even number of sites around each elementary plaquette show a dominant cpo = h/e 
periodicity. The lattices with an odd number of sites around each elementary plaquette show a 
dominant <^o/2 periodicity at vanishing diagonal disorder, with a positive weak localization- like 
magnetoconductance at infinite disorder fixed points. The magnetoconductance with both diagonal 
and off-diagonal disorder depends on the symmetry of the distribution of on-site energies. 
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I. INTRODUCTION 

The interplay between quantum fluctuations, corre- 
lations and quenched disorder is a notoriously difficult 
problem of solid state theory. Recently, however, a strong 
disorder renormalization group (RG) approach was in- 
troduced for a class a low-dimensional random quantum 
systems, for a review see[]J. The idea of this method 
has been introduced by Ma and Dasgupta[3 and has 
been widely used after Fisher had shownj^, Q| that the 
method provides asymptotically exact results for some 
one-dimensional problems. 

The applicability of the RG method is due to the na- 
ture of the low-energy fixed points in these systems. In 
the so called infinite disorder fixed points at larger scales 
the disorder grows without limits and therefore com- 
pletely dominates the low-energy, low-temperature prop- 
erties of the systems. In an infinite disorder fixed point 
both dynamical and static singularities are expected, de- 
scribed exactly by the RG method at low energy. The 
so called strong disorder fixed points, in which the disor- 
der has a limited strength, occur in different systems. At 
these fixed points only the dynamical singularities are 
expected to be calculated asymptotically exactly. Infi- 
nite disorder fixed points can be found in several ID sys- 
tems, for instance in the random transverse field Ising 
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chainQ, random antiferromagnetic Heisenberg champ), 
tight-binding model, equivalent in ID to the random XX 
model by a Jordan- Wigner transformation. Strong disor- 
der fixed points are characterized by Griffiths phases 0,0 
in the Heisenberg chain with randomly mixed ferro- and 
antiferromagnetic couplings 7] and in random quantum 
ladders 0. 

In higher dimensional systems we have less information 
about the properties of low-energy fixed points. Mod- 
els with discrete symmetry, such as the random quan- 
tum Ising model, generally have an infinite disorder fixed 
point |9j, whereas models with continuous symmetry, such 
as the random Heisenberg model have a strong disorder 
fixed point [Toj . For the random tight-binding model on 
the square lattice a logarithmically infinite disorder fixed 
point is conjectured [III H2- Generally, in higher dimen- 
sional systems the RG method can be implemented nu- 
merically and the analysis of the results is often difficult 
due to finite size effects. 

The core subject of our article is devoted to investi- 
gating the properties of the random tight-binding model 
with diagonal and off-diagonal disorder on fractal lat- 
tices through a strong disorder RG. We also make a 
comparison to random networks of conductors on Eu- 
clidean geometries treated by Green's functions, and to 
make contact with previous works on spin model, we dis- 
cuss also the random Heisenberg model on fractal lat- 
tices treated by a strong disorder RG. The Heisenberg 
and tight-binding models have different low-energy fixed 
points in ID and in 2D and we want to investigate the 
properties of this dimensional cross-over using different 
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types of fractal lattices. We use fractals with a sparse 
structure, so that the RG converges much more rapidly 
than on the compact higher dimensional Euclidean lat- 
tices. Besides the Sierpinski gasket for which we find 
strong disorder (i.e. not infinite disorder) fixed points, 
we study a family of quasi- ID lattices that can be cut 
into two pieces of arbitrary sizes by disconnecting only 
two bonds, for which we find infinite disorder fixed points. 
We also investigate the permanent current and the two- 
terminal conductance. We give a general argument show- 
ing that the permanent current is o /2-periodic for a 
tight-binding model with an odd number of sites around 
each elementary plaquette, with a symmetric diagonal 
disorder, and with a chemical potential in the middle of 
the band. We obtain an agreement between this gen- 
eral argument, the permanent current calculated from 
the RG, and the two-terminal current calculated from the 
RG and from Green's functions on Euclidean lattices. 

The structure of the paper is the following. The mod- 
els and the fractals are introduced in Sec[H]and the basic 
ingredients of the strong disorder RG method are pre- 
sented in Sec lIIII Our results are given in Sec II VI and 
IVl for the Heisenberg model and for the tight-binding 
model, respectively. We close our paper by a Conclusion. 
Some details are provided in the Appendices. 



(terms e k ): 

U = ^2 [t k j(4>)cf c k + £ i)fc (0)c£ci] - ^2e k c+c k . (3) 

(k,l) k 

The hopping integrals tk,i(<j>) are complex numbers and 
depend on the value of the magnetic flux <f> through an 
elementary plaquette: 

tk,l 0) = t k ,i exp J A.dv^j , (4) 

where the integral over the vector potential is from r k 
to r;, the two points corresponding to the lattice sites 
with labels k and I. The variables t k j correspond to the 
hopping integrals in the absence of a magnetic field. 

The on-site energies e k are chosen in the distribution 
V(e) = T a X*-> e max), where T is given by Eq. and 
the hopping variables t k> i are chosen in the distribution 
V{t) = T at (t, t max ). To avoid multiplying the number of 
parameters we restrict to the case a e — at and use the 
notation a = a e — a t . 



1. Correlation functions 



II. MODELS AND LATTICES 
A. Random Heisenberg model 

The random Heisenberg Hamiltonian is given by 

n = Y J Jfc,iS*.Sj , (i) 

(k,i) 

where the sum runs over all neighboring pairs of sites, 
and where S k corresponds to a spin-1/2. In the anti- 
ferromagnetic model we use < J k ^i < J ma x and the 
exchange distribution 

V{J)=T a {J,J^) = a^. (2) 

For the Heisenberg model with a symmetric distribution 
of exchanges we use — J m ax < Jk,l < Jmax- The absolute 
value of J is distributed according to Eq. @ . The moti- 
vation for introducing singular exchange distributions is 
that the strong disorder RG generates by itself singular 
renormalized coupling distributions even if we start from 
the regular distribution with a = 1. 



The correlation functions of the tight-binding models 
are defined by 



Ci(R) = M E^> ( 5 ) 

C2(R) = ttL£< c ^> 2 , (6) 

V '{k,i) 

with the constraint that the distance between sites k and 
I is equal to R. Here and in the following (...) stands 
for the thermal average (corresponding to the average in 
the ground state at T — 0), and the overline denotes 
averaging over quenched disorder. J\f(R) is the number 
of sites at distance R and A is a normalization constant. 

We first choose A in such a way that the histograms 
are normalized to unity: 

E^(i?) = L (7) 

The other possibility is to use unnormalizcd histograms 
corresponding to A = 1, that will be noted C[ 2 (-R)- C' 2 {R) 
will be used for comparing the conductance at different 
magnetic fields. 



B. Random tight-binding model 



2. Conductance, current and persistent current 



The random tight-binding Hamiltonian contains both 
off-diagonal disorder (terms t k j) and diagonal disorder 



The differential conductance Q(V, T, R) at a fixed tem- 
perature T, voltage V and distance R is defined as the 
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FIG. 1: Representation of the fractal lattices with n = 3, 
g = 4 (a); n = 4, g — 4 (b); n — 5, g = 5 (c); and the 
Sierpinski gasket with g — 6 (d) 



derivative of the two-terminal current X(V, T) with re- 
spect to the applied voltage: 



g{y,T,R) = —x(y,T). 



(8) 



For simplicity we carry out an averaging over all pairs 
of sites at distance R on the lattice, so that we cannot 
discuss the geometrical effects related to the position at 
which an extended contact has been connected to the 
fractal lattice. 

The persistent current </(</>, T) circulating in the sys- 
tem in the absence of an applied voltage is equal to the 
derivative of the energy U(cj>, T) with respect to the mag- 
netic flux: 



J{<t>,T) 



d_ 

80 



[/(«/», T). 



(9) 



C. The fractal lattices 

We use a family of fractal lattices based on regular 
polygons with n sides. We denote by g the number of 
generations. The lattices with n = 3, 4, 5 are shown on 
Fig. (a), (b) and (c). In the case n = 3 the number of 
sites N s and the number of links Ni is 



N s = 3 9 - 1 + - (10) 

N t - 3 9 . (11) 

In the case n > 3 we define A = (n 2 — An + 5)/[(n — 
l)(n - 3)], B = -n/[3(n - 3)], C = n/{n - 1), A' = 



(n-2)/(n-3). We have 
N„ = Ar 



3 a B + C 



Ni = A'n 9 + 3 9 B. 



(12) 
(13) 



The linear length is equal to 3 s so that the fractal dimen- 
sion is df — Inn/ In 3. 

We also use the standard Sierpinski gasket shown on 
Fig. ffl-(d). The number of sites and links are given by 



N s = 3 + - [39- 1 - 1] 



(14) 

n = o-. (15) 

The linear length is given by L = 2 9 so that the fractal 
dimension is df = In 3/ In 2. 

III. STRONG DISORDER RG APPROACH 

The strong disorder RG approach[l| as initiated by Ma 
and DasguptaQ works in the energy space and at each 
step a term in the Hamiltonian is eliminated to which the 
largest excitation energy (gap), denoted by fl, is associ- 
ated. During this decimation process the energy scale £1 
is lowered, and the number of sites of the lattice n(O), 
divided by the initial number of sites, is reduced. At 
the same time new couplings are generated between the 
remaining sites which were originally connected to the 
decimated sites. 

A. RG rules 

1. Heisenberg model 

During renormalization the Heisenberg model, if the 
lattice is more complex than the linear chain, it will 
transform into a set of spins of different size between 
them if there are both ferromagnetic and antiferromag- 
netic couplingsjll. The renormalization rules, as given in 
the literature 0, 0, Q] involve two processes: spin- 
cluster formation and singlet formation. Two spins S\ 
and 52 coupled by a strong exchange | Ji,%\ ~ give rise 
to a spin S = Si + S2 if Ji,2 < (ferromagnetic coupling) 
and to a spin S = |5i — S2I if Ji,2 > (antiferromag- 
netic coupling) and Si 7^ S2. During this spin-cluster 
formation process the renormalized couplings are calcu- 
lated by first order perturbation theory. There is singlet 
formation if Si = S2 are coupled antiferromagnetically, 
in which case the pair of spin Si and S2 is eliminated and 
the renormalized couplings are calculated by second order 
perturbation theory. The expression of the renormalized 
couplings can be found in the literature fH lei ITolll4| . 

2. Tight-binding model 

In the case of the tight-binding model there are two 
transformations, (i) If the largest gap of the system is 
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associated to a diagonal term, say = Q, a fermion or a 
hole is frozen at site k, and this site is decimated out. We 
will refer to this process as an e-transformation. (ii) If the 
largest gap of the system is associated to an off-diagonal 
term, say tk.i = 0, a fermion is frozen in a "dimer" made 
of sites k and I and this dimer is decimated out. We 
will refer to this process as a ^-transformation. The RG 
transformations for a bipartite lattice with = 0, for 
all k involves only the ^-transformation and can be found 
already in the literature^^- A specificity of the non 
bipartite lattice RG transformations given in AppendixlAl 
is that diagonal disorder can be self-generated by the RG, 
even though diagonal disorder is not present in the initial 
condition. 



B. RG flow 

During renormalization we monitor the distribution 
function of the different parameters of the Hamiltoni- 
ans in Eqs. and as a function of the energy scale 
n. In different systems the distribution of couplings, say 
P(J, might have different characteristics as the fixed 
point, fl* = 0, is approached. 



1. Infinite disorder fixed point 

At the so called infinite disorder fixed points the dis- 
tributions broadens without limits, so that the ratio of 
typical renormalized couplings goes to zero or to infinity. 
In this fixed point it is convenient to introduce the log- 
energy scale, T = — ln(no/f2) (where fio is a reference 
energy), related to the length-scale, L ~ ^(r)^ 1 ^, as 



(16) 



where n(Y) is the density of non-decimated points. The 
exponent if) is one of the exponents characterizing the 
universality class of the infinite disorder fixed point. 

This type of fixed point is realized, among others, for 
the random antiferromagnetic Heisenberg spin chain and 
for the random tight-binding model in ID, for which %j) = 
1/2 in the two cases. In an infinite disorder fixed point 
the distribution of the lowest gap A follows the functional 
form: 



P L (A)dA = P 



In A 



L'-' 



dlnA 



(17) 



2. Logarithmically infinite disorder fixed point 

At this fixed point the relation between log-energy 
scale and length-scale is: 



(InL) 



UJ+l 



(18) 



thus the exponent ip is formally zero. This type of be- 
havior is found in the random tight-binding model on the 



square lattice. Gade[TT| found w = 1, whereas Motrunich 
et al.[v^ conjectured u = 1/2, see also Ref . [l3| . At this 
fixed point the distribution of the lowest gap is given by: 



P L (A)dA = (lnL) w+1 P 



In A 



L(lnL)- 



dlnA 



(19) 



3. Strong disorder fixed point 

At this fixed point the relation between energy scale 
and length-scale is: 



(20) 



where z is the disorder induced dynamical exponent. To 
obtain the true dynamical exponent zq of the system, 
one should also consider deterministic quantum fluctu- 
ations which could results in a vanishing gap and thus 
a quantum dynamical exponent z q . In the non-random 
system z — z q , whereas in the presence of quenched dis- 
order zq = max(z 9 ,z). Thus, for z > z q the dynamical 
singularities of the system are determined by disorder 
effects, otherwise disorder will results in confluent (sub- 
dominant) singularities. (In the infinite and logarithmi- 
cally infinite disorder fixed points z is formally infinite.) 

Strong disorder fixed points are found, among others 
in random antiferromagnetic Heisenberg models on the 
square and cubic lattices ^(J and generally this behavior 
is related to quantum Griffiths phases0,|a,|(|. 

At this fixed point the distribution of the lowest gap is 
given by: 



P L (A)dA = L- Z P 



A 
I 7 



dA 



(21) 



If the excitations are localized, then for a small, fixed A, 
Pl(A) is proportional to the volume, i.e. Pr,(A) ~ L d , 
from which follows in the limit large L limit: 



P(A) 



-l+d/z 



(22) 



This relation is used in numerical calculations to de- 
termine the dynamical exponent. The asymptotic distri- 
bution of the gaps at a strong disorder fixed point has a 
power-law form, that should be compared to the initial 
distribution in Eq.Q. There are random quantum sys- 
tems in which z depends on the initial disorder and z is 
approximately proportional to I /a. This can be found, 
among others for random Heisenberg spin ladders 0- In 
other systems with frustration, z is disorder indepen- 
dent (see the three-dimensional spin glass models Jll§), 
and large spins are formed in these cases, as discussed in 
the following section. 



4- Singlet and large spin formation 

In random Heisenberg models during decimation the 
typical size of the generated new spins can be of two 
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types, i) The renormalized spins are typically 5=1/2 
(or O(l)) and do not grow under renormalization. In this 
case in the ground state typically singlets are formed and 
we refer to this process as random singlet formation. This 
type of behavior is found in the ID random antifcrromag- 
netic Heisenberg model 0. ii) Another possibility is that 
during renormalization large spins are formed, the size of 
which continuously grows during renormalization. The 
typical size of the effective spin, S e ff, is asymptotically 
related to the number of decimated spins, N as: 



S, 



eff 



(23) 



The characteristic exponent £ would be C = 1/2 in an 
uncorrelated decimation approximation 0, ll(1 |. Large 
spin formation can be found, among others in the 
±J Heisenberg models in one- |3, two-HU, and three- 
dimensions 



C. RG calculation of correlations and transport in 
the random tight binding model 

1. Correlation functions 

During renormalization, as f2 is gradually decreased, 
some electrons and holes becomes localized at given sites 
and do not contribute to the correlation functions in 
Eq.JSJ) and ©. There is, however, a complementary pro- 
cess by which electrons arc bound on delocalized orbitals 
connecting remote sites, say k and I. Therefore correla- 
tions between k and / are 1/2, if tki at the decimation 
energy, \tki\ — fi, is positive and it is —1/2 otherwise. 
These rare regions dominate the average correlation func- 
tion. Typical correlations, measured between two sites 
that are not involved in a delocalized bound state, are 
much weaker and scale as ~ tki/Q- This is indeed negli- 
gible if the distribution of tjy is broad. Thus we obtain for 
the leading behavior of the average correlation function, 



Ci(to,R) 



X 



2N(R) 



E 

(k,i) 



Ski 



(24) 



where the sum is over all frozen pairs at distance i?, and 
Ski = sign(ife(). We have a similar expression also for 

c 2 (n,R). 



same rare regions are dominant for the average unnor- 
malized correlation function C' 2 {J^ , R) we obtain the rela- 
tion: 



g(V,T,R) = -C 2 (Sl(V,T),R) 



(25) 



A detailed derivation of Ea. l(^5|l . based on the perturba- 
tive theory of transport processes can be found in Ap- 
pendix [5] 

The total current J to t ( V, T, R) in Eq. JSJ is obtained by 
integrating the differential conductance (|25J) over voltage, 
which for a small fixed temperature, T <C eV , is equiva- 
lent to an integration over the energy scale, f2: 



X to t(V,T, R)k 



fi(V,T) 



g{n',R)dn' 



(26) 



Eq. I|26|) constitutes a generalization of the Kubo formula 
since Xt t(V, T, R) is not necessarily linear in V in the 
limit V — > 0. The differential conductance for a large sep- 
aration between the contacts corresponds to low energies, 
so that we use the standard form of the RG asymptoti- 
cally exact in the low energy limit at an infinite disorder 
fixed point. 



3. Persistent current 

The zero temperature persistent current in Eq.@ is 
related to the average energy 



U(</>,T) 



n 

fl(T) 



(27) 



where p*(f2) and p e (f2) are the probabilities to perform 
respectively a t- or and e-transformation at energy scale 
fi. We have subtracted in Eq. I|27l) the energy U(4>, +oo) 
at infinite temperature. The infinite temperature energy 
in a t transformation is vanishingly small since the spec- 
trum of a dimer is symmetric. The infinite temperature 
energy in an e-transformation is equal to e/2. The e state 
is occupied only below the Fermi level, leading to the en- 
ergy c6(-c) (with 0(e) = 1 if e > and 0(e) = if e < 0). 
Subtracting e/2 yields -|e|/2, leading to Eq. (f?7|) . 



IV. RANDOM HEISENBERG MODELS 



2. Conductance and current 

Here we consider the system with a small applied volt- 
age V and at a small temperature T. This will generate 
an energy scale in the system, Q(V,T) — max(eV r , T), 
and the renormalization process stops at £1 = Q(V,T). 
The average value of the differential conductance is dom- 
inated by rare regions, in which, at energy scale Q(V,T) 
a frozen extended electron state connects the two elec- 
trodes, that are separated by a distance R. Since the 



To make contact with previous works on localized spin 
models in ID, 2D and 3D we first consider the Heisenberg 
model on fractal lattices before investigating the tight- 
binding models on the same lattices. In the numerical 
calculations we considered random Heisenberg models on 
the family of fractals with n = 3, ...10 and on the gasket, 
for a purely antiferromagnetic exchange distribution and 
for the symmetric exchange distribution with mixed fer- 
romagnetic and antiferromagnetic couplings. The strong 
disorder RG procedure is implemented numerically and 
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FIG. 2: (a) Collapse plot of the gap spectra for the antifer- 
romagnetic Heisenberg model on the fractal with n — 3. (b) 
Collapse plot of the gap spectra for the Heisenberg model with 
a symmetric distribution of exchanges (noted ±J models on 
the figure) on the fractal with n = 3. We use a — 0.25 (□, 
red) and a = 0.5 (•, green) on (a) and (b) respectively. We 
collapsed the gap spectra with g — 3,4, 5, 6, 7 generations, (a) 
corresponds to the scaling 117H with ip — 0.44 while (b) corre- 
sponds to the conventional disorder scaling given by Eq. 1211 1. 
The corresponding values of z are shown in Fig. [3] The gap 
spectra of the antiferromagnetic Heisenberg model and of the 
Heisenberg model with a symmetric distribution of exchanges 
are shown on (c) and (d) for g = 3 (■, red), g — 4 (□, green), 
g = 5 (o, blue), g = 6 (•, purple), and g = 7 (A, blue). 



the decimation is performed either up to the last re- 
maining (large) spin or up to the last singlet and the 
corresponding gap, A, is calculated. From the distribu- 
tion of the gap -P(A), the properties of the fixed point 
are determined using the results overviewed in Sec lIIIBl 
The number of generations is g — 3, 4, 5, 6, 7 for n = 3, 
g = 3, 4, 5, 6 for n = 4, g = 3, 4, 5 for n = 5, 6, g = 3, 4 for 
n = 7, 8, 9, 10. We use a = 0.25, 0.5, 0.75, 1 for n = 3, 4, 
and a = 0.25 for n = 5,6, 7. We use 250000 realizations 
of disorder in each simulation. 



A. Distribution of the gap 

1. Known results on regular lattices 

In ID the random antiferromagnetic Heisenberg 
model .4] has an infinite disorder fixed point with the ex- 
ponent i\) — 1/2. On the contrary the ±J model in ID 
has a strong disorder fixed point Q with a finite dynam- 
ical exponent, z — 2.3 and a large-spin with £ = 1/2. 
Quasi-one-dimensional models, i.e. spin ladders with 
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FIG. 3: (a) ip versus n for the antiferromagnetic Heisenberg 
model, (b) z versus 1/q (proportional to the strength of disor- 
der) for the Heisenberg model with a symmetric distribution 
of exchanges (noted ±J model on the figure). The curve (1) 
(□, green) corresponds to n = 3 and n — 5 with z determined 
from the scaling form 121H of the gap spectra (the points for 
n = 3 and n — 5 are superimposed). The error-bars are com- 
parable to the size of the symbols. The curves (2) (■, blue) 
and (3) (•, purple) correspond to n — 3 and n — 5 respec- 
tively, with z determined from the slope of the rare event tail 
of the gap spectra, (c) Variations of z (□, red) and z/df 
(■, green) versus n for Heisenberg model with a symmetric 
distribution of exchanges with a = 0.25. (d) Variation of z 
(determined from the slope of the rare event tail of the gap 
spectra) versus g for the antiferromagnetic Heisenberg model 
(□, red) and the Heisenberg model with a symmetric distri- 
bution of exchanges (■, green) on the gasket with a = 1. 



random antiferromagnetic couplings have a strong dis- 
order fixed point in which the dynamical exponent is ap- 
proximately proportional to the strength of disorder, i.e. 
to a -1 starting with the distribution in Eq.Q 8]. 

In 2D the random antiferromagnetic model on the 
square lattice, which is non-frustrated, has a strong dis- 
order fixed point, in which the dynamical exponent is 
approximately constant for weak disorder [lfj, but for 
stronger disorder z becomes disorder dependent 0]. On 
the other hand for frustrated models (of geometric ori- 
gin or for the ± J model) there is a strong disorder fixed 
point with a large spin. This fixed point with z = d and 
C = 1/2 is disorder independent and was called a "spin 
glass" fixed point 01 . There are less accurate numerical 
results in 3D, but it is likely that a similar scenario as in 
2D holds also in 3P[!o|. 
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2. Random antiferromagnetic model 

The results shown in Figs. |21 and |21 suggest that for 
the antiferromagnetic Heisenberg model the fractals with 
n even correspond to a infinite disorder fixed point de- 
scribed by the usual scaling with ip = 1/2 whereas the 
fractals with n odd correspond to non trivial infinite dis- 
order fixed points with ip lower than 1/2 (see Fig.|3J-(a)). 
The ground state is always a singlet for n even whereas 
there is large spin formation for n odd. An inspection of 
the RG flow indicates that spins 1 and higher are gen- 
erated in the RG for n even but correspond to an irrel- 
evant perturbation since they disappear soon after they 
are formed. The uncertainty in the determination of ip 
for n = 7,8,9, 10 on Fig. |3] is due to finite size effects 
since we cannot simulate large values of g for these val- 
ues of n. The error-bar is estimated from the finite size 
effects with n = 4. 



+J model, n=3, £=0.43 
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3. Random ±J model 

The gap spectra of the Heisenberg model with a sym- 
metric distribution of exchanges are well described by the 
strong disorder scaling J2J, for both n even and n odd. 
The value of z can be extracted either from the collapse 
of the gap spectra given by Eq. I|21ll or from Ea. l|22[) by 
fitting the rare event tail. We have performed both anal- 
ysis and obtained - within the error of the calculation - 
identical results. The disorder induced dynamical expo- 
nent is found proportional to the initial strength of disor- 
der 1/a, and we obtained approximately: z ~ 0.3df/a. 
This corresponds to the green dashed line in Fig. O (c) 
for a = 0.25. 



FIG. 4: Collapse plot of the distribution of the spin in the 
ground state for the Heisenberg model on the fractal with 
n — 3, for g = 4, 5, 6, 7 and £ = 0.43. The collapse plot 
corresponds to Eq. 1231 . 



scaling in Eq. (|23|l is well obeyed. We find for the expo- 
nent £ = 0.43 ±0.02 for the antiferromagnetic Heisenberg 
model with n — 3 (see Fig. and £ = 0.5 in all other 
cases. 



Conclusion 



4- Results on the Sierpinski gasket 

For the Sierpinski gasket we find that disorder at the 
end of the RG decreases as the system size increases, 
both for the antiferromagnetic Heisenberg model and the 
Heisenberg model with a symmetric distribution of ex- 
changes (see Fig. E2-(d)). The effective, size dependent 
dynamical exponent seems to tend to zero as 1/lniV. 
This implies that either i) there is a small finite gap in the 
system or ii) although the gap is vanishing the dynami- 
cal correlations decay faster than a power-law. However, 
due to the limited sizes of the systems we cannot further 
study these scenarios. 



B. Distribution of the spin of the ground state 

There is large spin formation for the antiferromagnetic 
Heisenberg model on the fractal with n odd, and for the 
Heisenberg model with a symmetric distribution of ex- 
changes for all values of n. In all cases we find that the 



The empirically obtained critical behavior of the ran- 
dom Heisenberg model on different fractal lattices is sum- 
marized in Table J] The low-energy fixed point seems to 
depend on the form (symmetry) of disorder, the symme- 
try of the fractal lattice, but there is no direct relation 
to the fractal dimension of the lattice. Interestingly the 
studied fractal lattices constitute low-energy universality 
classes that can not be found in the regular ID and 2D 
lattices. 



V. RANDOM TIGHT-BINDING MODELS 

The detail of the RG transformations of the tight- 
binding model can be found in Appendix The 
interferences between different trajectories in the t- 
transformation arc eliminated if the sum of the different 
terms in the renormalized coupling is replaced by their 
maximum, an approximation used in the study of spin 
models. To keep interferences, we do not use the maxi- 
mum in the expression of the renormalized couplings. 
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TABLE I: Low-energy fixed points of random Heisenberg models on fractals. For a comparison results on regular ID and 
2D (square) lattices are also presented. (AF: antiferromagnetic couplings, ±J: mixed antiferromagnetic and ferromagnetic 
couplings, ID: infinite disorder fixed point, SD: strong disorder fixed point, RS: random singlet formation, LS: large spin 
formation.) (+) For n = 3, C < 1/2. 





n — odd 


n = even 


gasket 


ID 


2D 


AF 


IDty = 1/2) 
LS(C = l/2 (+) ) 


lD(tp < 1/2) 
RS 


SD(*-0)(?) 
LS(C = 1/2) 


ID(^ = 1/2) 
RS 


SD 
RS 


±J 


SD(z ~ 1/a) 
LS(C=l/2) 


SD(z ~ 1/a) 
LS(C = 1/2) 


SD(*-0)(?) 
LS(C = 1/2) 


SD(z = 2.3) 
LS(C = 1/2) 


SD(z = d) 
LS(C = 1/2) 
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FIG. 5: Fraction of frozen sites n E (defined by Eq. 1281 1 versus 
/ (a) for the fractal with n = 3 and g = 4 ■ red), g = 5 
(•, green) and g = 6 (A, blue); and (b) for the gasket with 
with g = 4 (■, red), <? = 5 (•, green) and g = 6 (A. blue) 
generations. We used a = 1. 



A. Features of the RG flow without extrinsic 
diagonal disorder 

1. Fraction of frozen sites 

We first consider the case where diagonal disorder is 
not introduced in the initial condition. We thus start 
with efc = for all k. For a bipartite lattice (for in- 
stance the fractal with n even) no diagonal term is gen- 
erated during the RG transformations, which can be seen 
from the decimation rules given in Appendix lAl On the 
contrary for a non bipartite lattice (for instance for the 
fractal with n odd) also non vanishing values of ejt are 
generated, and at the end of the RG an average number 
N e of the total number of sites N s has been frozen by 
the e-transformations. We have calculated the fraction 
of frozen fermions. 



NJN S 



(28) 



for both the fractal with n = 3 and for the gasket. We 
work in the grand canonical ensemble where the chemical 
potential is fixed. 

We show in Fig. [3] the variation of n e as a function 
of / = (f>/4>Q- First we note that n t decreases slowly 
with increasing system size: in the limit N s — > oo our 
data indicate an approximately logarithmic dependence, 
n e ~ l/hxN s . In a finite size system n € shows <fco/2- 
oscillations so that the number of frozen fermions van- 
ishes at / = 0.25,0.75,1.25,... and there is a tendency 
of the system to be more localized for / = 0,1/2,1..., 
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FIG. 6: Average coordinance of the graph of "active" spins 
as a function of F, for different values of Fo, for the fractal 
with n = 3. (a) corresponds to g = 6 generations and / = 0. 
(b) corresponds to g = 7 and / = 0. (c) corresponds to g = 6 
and / = 0.25. (d) corresponds to g — 7 and / = 0.25. The 
different curves in each panel correspond to To = 10 (■), 
Fo = 30 (•) and Fo = 50 (a). The red color is used for g = 6 
and the green color for g — 7. We used a = 0.25. 



in agreement with the positive magnetoconductance dis- 
cussed later in section IV El The cancellation of n £ for 
/ = 0.25, 0.75, 1.25, ... can be understood from the renor- 
malization of e in the f-transformation (see Eq. IA3|) . 
If the sites "1" and "2" are eliminated in a three-site 
cluster {1,2,3} then the interference of the trajectories 
3 — > 1 — > 2 — > 3 and the trajectory 3 — > 2 — > 1 — > 3 are 
such that £3 = 0, so that no diagonal disorder is gener- 
ated by the RG for these values of / = 
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FIG. 7: Variation of the average clustering coefficient c'^ v (T), 
see text, as a function of V, for different values of To, for the 
Euclidian 2D square lattice with / = 0. (a) corresponds to 
L = 16 and (b) corresponds to L — 28. The different curves 
in each panel correspond to To = 10 (■, red), To = 30 (•, 
green) and Fo = 50 (A, blue). We used a — 0.25. 



2. Average coordinance 

A useful information is encoded in the connectivity of 
the graph of "active" sites (the sites that have not been 
eliminated at the energy scale f2). More precisely we con- 
sider the connectivity of clusters of "strongly connected" 
sites, which are defined by having hopping amplitudes 
tk,i that exceed a limiting value: 



In 



<r , 



(29) 



where To is a given cut-off. 

We calculate the average coordinance c av of this graph 
as a function of the log-energy scale T, where c av is de- 
fined as the ratio of the average number of bonds and 
the average number of sites. The average coordinance 
c av is expected to start from a value of order unity, first 
increase, reach a maximum, and then decrease to zero at 
large T since only a few sites are surviving at the end of 
the RG. 

Cav(F) is shown in Fig. for the fractal with n = 3, 
for two sizes corresponding to g — 6,7, and for / = 
and / = 0.25. It is visible that there are practically no 
finite size effects for g — 6, 7. The average coordinance is 
larger for / = 0.25 than for / = 0, compatible with the 
existence of a negative magnetoresistance discussed in 
section IV El the system is more connected and therefore 
more delocalized for / = 0.25. 

By comparison we carried out a similar simulation for 
the Euclidian two-dimensional square lattice. The finite 
size effects indicate that after a transient the graph of 
surviving sites becomes extremely connected whereas it 
is loosely connected in the case of the fractal with n = 3. 
This geometrical effect is related to the fact that the al- 
gorithm becomes extremely slow in the case of the Eu- 
clidian square lattice as compared to the fractal lattices. 
More quantitatively it is convenient to evaluate the clus- 
tering coefficient c' av (T), which is the average fraction 



10" 



10"' 



10"' 



o 
_i 

ST 



10" 



10"' 



10"' 



10"' 



n=3, f=0 

- to f % (a) - 

J A 

t i a 

- $ • A 

III ' A 

u \ \ 


« n=3, f=0.25 
*\ < b > " 

T& * \ 

\ ' * V 
l\ \ 

— i — i h r 


— h- — I 1 — 

h u n=4, f=0 

*A**V% (c) - 

i 5 i • N 

1 • ^ 

i i i 


*v n=4, f=0.25 

-m % (d) - 
■if* \ 

If" • A 
" I"? » 'a 

i A 
• 

1 1 1 



10 20 30 10 20 30 40 
-Log(A) 

FIG. 8: Gap spectra of the tight binding model with (n — 3, 
/ = 0) (a), (n = 3, / = 0.25) (b), (n = 4, / = 0) (c), and 
(n = 4, / = 0.25) (d). The simulations on (a) and (b) with 
n = 3 correspond to for g = 4 (■, red), g = 5 (o, green) g = 6 
(•, blue), and g ~ 7 (A, purple). The simulations on (c) and 
(d) with n — 4 correspond to g = 3 (■, red), g = 4 (o, green), 
3 = 5 (•, blue), and g = 6 (A, purple). We used a = 1. 



of bonds, relative to the maximum number of possible 
bonds, N S (N S — l)/2. For sufficiently large values of Fo 
the quantity c^ v (r) reaches a plateau at around 0.7, see 
Fig0 indicating that the graph of active sites becomes 
strongly connected. 



3. Gap spectra 

We show in Fig.[S]for the fractals with n = 3,4 the gap 
spectra of the random tight-binding model with e max = 
and / = 0, 0.25. In all these cases there is infinite disorder 
scaling 16], see Sec lIIIB ll The collapse plots according 
to Eq. il7|) show that tp = 0.54 ±0.04 for (/ = 0, n = 3), 
ip = 0.5 ± 0.03 for (/ = 0, n = 4), if) = 0.48 ± 0.02 for 
(/ = 0.25, n = 4). To obtain the scaling for (/ = 0.25, 
n = 3) we repeated the simulations in Fig.[SJ-(b) without 
the last e-transformation and found that the scaling is 
well described by = 0.48 ± 0.02. 

The gap spectra for the gasket shown in Fig.[!5]indicate 
strong disorder scaling, see Sec llll B"3l The parity effect 
for / = 0.25 is similar to the fractal with n = 3 and 
/ = 0.25 in Fig. |SJ see0|. For the Sierpinski gasket the 
number of sites N s is odd for g odd and even for g even. 

We show in Fig. ^| the variation of the ratio z/df 
versus / for the Sierpinski gasket, as calculated from 
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FIG. 9: Gap spectra of the tight binding model on the Sier- 
pinski gasket with / = (a) and / = 0.25 (b) and with g = 3 
(red, □), g = 4 (green, ■), g = 5 (blue, o), g = 6 (purple, •) 
and g = 7 (blue, A) generations. We used a — 1. 
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FIG. 11: Normalized correlation function C2(T,R) versus R 
for the fractal with n = 3 with g — 7 generations, a — 0.25 
and / = (a), / = 0.25 (b). The different curves correspond 
to r = 5 (■, red), F = 30 (•, green). The solid lines corre- 
spond to the fit igOJ, with x = 1.6 and 5(F) = 2.5 for T = 5, 
and £(r) = +oo for T = 30. 
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FIG. 10: Ratio a/d/ versus / for the Sierpinski gasket with 
(g = 4, a = 1) (■, red), (<? = 5, a = 1) (•, green) and (<? = 6, 
a = 1) (A, blue), (g = 4,q = 0.25) (□, red), (3 = 5,a = 0.25) 
(o, green) and (g = 6, a = 0.25) (A, blue). The ratio z/df is 
smaller than unity for sufficiently large system sizes. 



Eq. (|22p. For a weaker initial disorder (a = 1) the dy- 
namical exponent is practically independent of / and has 
no finite size effect, so we can estimate z/df — 0.5. The 
finite size effects are also stronger for stronger initial dis- 
order (a = 0.25). The dynamical exponent for very large 
system sizes is probably /-independent, and the appar- 
ent variations are within the accuracy of the numerical 
calculation. From the available data it is difficult to es- 
timate the limiting value. However, z/df (the strength 
of disorder at the fixed point) increases with 1/a (the 
strength of disorder in the initial condition), so that we 
might assume that the limiting value of z/df could be 
close to 0.5 for a = 0.25. 



The definition of the correlation function is given in 
section lill C II The correlation functions obtained from 
the simulations are well fitted by 



1 



C 2(r,J?)~^exp 



R 



(30) 



with a correlation length £(r) diverging like £(T) ~ r 1 ^ 
at an infinite disorder fixed point, see Ea. Hlbj) . For the 
Euclidean ID random tight-binding model we have x = 1 
from the simulation, in agreement with the analytical 
prediction Q. For the fractals with n = 3 and n = 4, we 
obtain the same values of the decay exponent x and the 
correlation length £(T), for Ci(r,i?) and C 2 (T,R). This 
is expected since the dominant contribution to the cor- 
relation functions is due to the same type of rare regions 
for both functions. The decay length of Ci(T, R) is equal 
to the decay length of C2(F, R), a feature of ID Euclidean 
systems that we recover on the sparse fractal lattices. In 
usual disordered systems with diagonal disorder and no 
off-diagonal disorder, the decay length of C\(R) at zero 
temperature is equal to the elastic mean free path, and 
the decay length of C2(R) at zero temperature is equal to 
the localization length. At finite temperature there exists 
an additional exponential decay over the thermal length 
£r- The situation considered here with off-diagonal dis- 
order and diagonal disorder can be viewed as a tempera- 
ture dependent mean free path equal to the temperature 
dependent localization length. Like in ID systems, there 
is no genuine diffusive regime in which the localization 
length would be much larger than the elastic mean free 
path. Interestingly, we still obtain a weak localization- 
like magnetoresistance for the fractal with n = 3 (see 
Sees ESI and EH . 

The variation and the fit of C2(I\ R) versus R is shown 
in Fig. ^] for the fractal with n — 3. We used x = 1.6 
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FIG. 12: Variation of the fraction of sites n e frozen by an 
e-transformation as a function of e max , the upper cut-off in 
the distribution of e, for the ID chain (a) and the fractal 
with n = 4 (b). For the ID chain on (a) we used chains of 
length L = 100 (■ and □, red), L = 400 (• and o, green) 
and L = 1600 (A and A, blue). For the fractal with n — 4 
on (b) we used g — 3 (■ and □, red), g = 4 (• and o, green) 
and g = 5 (A and A, blue). The open symbols correspond to 
a = 1 and the filled symbols correspond to a — 0.25. 

for all values of a but, as expected, £(r) increases as a 
increases (not shown in Fig. Illjl . The same value of x is 
obtained for n — 4 (not shown in Fig. Illfl . We find no 
significative difference in the values of x and £(r) between 
the cases / = and / = 0.25. 

There exists a parity effect like in ID. Namely the 
correlations between two sites at distance R are much 
stronger if R is odd. The origin of this effect in ID is 
that all the sites in between should be coupled in local- 
ized pairs in order to have an extended localized orbital 
between the two remote sites at distance R. This is ob- 
viously only possible for R odd. Similar effect might be 
due to the ID chains embedded in the fractal with n = 3. 



C. Effect of extrinsic diagonal disorder at / = 

1. Fraction of frozen sites 

We now consider the case of a finite extrinsic diago- 
nal disorder with / = 0. Increasing e max tends to in- 
crease the number of e-transformation and thus the frac- 
tion of frozen sites 7"i e (e max ). For bipartite lattices, such 
as the ID chain and the fractal with n — 3 the lim- 
iting value of n. e (e max ) is finite in the limit e max — > 0, 
even for a finite lattice. This behavior is illustrated in 
Fig. El The limiting behavior of n e (e max ) can be esti- 
mated as follows. For a small fixed e max , the RG trans- 
formation goes along the trajectory of the e max = case, 
until the log-energy scale T = ln(eo/e max ) is reached, 
with some reference value eo- At this point n e (e max , N s ) 
is proportional to the number of sites N s . The infi- 
nite disorder scaling given by Eq. I|16f) corresponds to 
ln(eo/e max ) ~ T ~ Nf^ df . We thus obtain a logarith- 
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FIG. 13: Variation of the fraction of sites n e frozen by an e- 
transformation as a function of e m ax, the upper cut-off in the 
distribution of e, for the fractal with n = 3. We used g = 4 (• 
and o, green) and g = 5 (A and A, blue). The open symbols 
correspond to a = 1 and the filled symbols correspond to 
a = 0.25. We obtained a similar behavior for n — 5. 

mic dependence as n e (e max ) ~ (ln(e /e max )) <i ^ ,/; . The 
numerical data in Fig. ll2l are compatible with this form, 
although a simple power-law form n e (e max ) ~ e^ax' with 
an effective exponent y = 0.14 seems to work, too, both 
for the ID chain and the fractal with n = 4. 

For a non-bipartite lattice, such as the fractal with 
n = 3, n e tends to a constant in the limit e max — > for a 
finite lattice, see Sec IV A ll This is illustrated in Fig. 1131 
For a small e max , n e increases as the system size increases 
for bipartite (the ID and the fractal with n = 4) lattices 
and decreases for non-bipartite (the fractal with n = 3) 
lattices. 



2. Correlation length 

The average correlation length £ a u(r, e max ) is calcu- 
lated from the average correlation function C2(T, R). Us- 
ing the functional form in Eq. I|30f) we fix the decay ex- 
ponent x to its value for e max = (see Sec. IV B|) . and 
then extract the value of £ av (T, e max ) from the numeri- 
cal data. The numerical results are shown in Fig ll4l for 
the ID chain (a) and for the fractals with n = 3 (b) and 
n = 4 (c). In the limit T ^> ln(e roax /eo) the average cor- 
relation length has a T independent limiting value which 
depends on ln(e max /eo) and the functional form seems to 
be linear for all type of lattices. Here we can observe 
an analogy with the behavior of the average correlation 
length in a disordered Heisenberg chain used as a model 
of disordered spin-Peierls system at a finite temperature 
T, as calculated in Ref. 18]. Making the correspondence 
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FIG. 14: Average correlation length £au(F,e max ) as a function of e ma x for several values of F, for the ID chain (a), the fractal 
with n = 3 (b) and the fractal with n = 4 (c). In (a) we use r = 7 (■, red), r = 9 (•, green), F = 11 (A. blue). In (b) and (c) 
we use F = 5 (■, red), F = 8 (•, green), V = 11 (A. blue), We use a = 0.25, L = 6400 for the ID chain, g = 7 for the fractal 
with n = 3, and g = 6 for the fractal with n — 4. Fhe fits correspond to £av(r,e max ) = a In (e max /f2o) + b, with a = 0.35 and 
b = 1.2 for the ID chain, and a = 1 and 6 = 1.8 for the fractals with n = 3 and n = 4. We obtained qualitatively similar results 
for other values of a, and for the "±e" model with a symmetric distribution of e. 



with the energy scales in the two problems, e max «-> T, 
we arrive to the relation £ av (F, e max ) ~ ln(e max /e ) for 
F 3> ln(e max /eo), as observed numerically. 

D. Persistent current 

1. General argument 

We expect the averaged persistent current to exhibit a 
</>o/2 periodicity for any tight-binding lattice with an odd 
number of sites around its elementary plaquettes, pro- 
vided the distribution of diagonal disorder is even and the 
chemical potential fi is fixed equal to zero. The follow- 
ing argument is inspired by the discussion given in 19] for 
Aharonov-Bohm current oscillations in a mesoscopic ring. 
The eigenvalue equation for the single-particle Hamilto- 
nian © reads: 

Eip k = -e k ipk + E tkl ( 31 ) 
i 

If we multiply this equation by —1, we get an eigenstate 
with energy —E, for a tight-binding model where on-site 
potentials and hopping amplitudes have been turned into 
their opposite values. On a bipartite lattice, changing 
tki{4>) into —tki(4>) is simply equivalent to performing a 
gauge transformation, which leaves the external magnetic 
field unchanged. But on a lattice (such as a triangular 
lattice) with elementary odd cycles, this change of sign 
in the hopping amplitudes is equivalent to changing the 
elementary flux <f> into <p + 4>o/2. Denoting single-particle 
energy eigenvalues by E a (<p, {efc}), we have then 

E Ea(4>, {e k })6{n-E a ) (32) 

a. 

+ ]T E {cj> + ^-,{-e k })d(E p + n) = (33) 



Let us then average these expressions over disorder real- 
izations. Since the e k distribution is assumed to be even, 
this yields: 



E E a {4>)9^ -E a ) + J2 M4> + y )^/3 + m) = o 

a 

The sum-rule on the single-particle spectrum yields 



E W + Y mE ? + M) + e{ -^ - E ^ 
= -E^=° ( 35 ) 

k 

So finally: 



E E <* (MO* - E «) = E W + y W-a* - E e) 

a 3 

(36) 

For fi = 0, this implies the </>o/2 periodicity of the sys- 
tem's total energy. The above argument also holds for 
the fractal considered here with n — 3 or for the Sierpin- 
ski gasket even though they have also some elementary 
loops with an even number of sites. This is because these 
larger loops enclose an area which is an even multiple of 
the elementary triangle area. So changing <fi into 4>+4>o/2 
does not change the total flux modulo <po in larger ele- 
mentary loops, which is compatible with a sign reversal 
of all hopping amplitudes. 

A direct experimental test of these predictions is how- 
ever difficult, since we have considered single channel 
tight-binding models, and it is not possible to put a flux 
quantum through an area of atomic size. The qualitative 
difference between bipartite (n = 4) and non-bipartite 
(n = 3) lattices also manifests itself in conductance oscil- 
lations (see sections IV El and IV Ffl . As shown in Ref. 0, 
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FIG. 15: Evolution of the harmonics Ai(e max ) (open symbols,red) and ^(emax) (filled symbols, green) of the permanent 
current, as a function of e max . Panel (a), (b) and (c) correspond respectively to the fractals with n = 3, n — 4 and to the 
Sierpinski gasket. The different curves within each panel correspond to the "±e" model with a symmetric distribution of et (□ 
and ■), to the model with et > (o and •), and to the "— e" model with et < (A and A > 



our argument for 0o/2 periodicity in lattices with odd 
elementary cycles also applies for the conductance calcu- 
lated using the Kubo formula. 

2. RG results 



tices regarding the behavior of the first harmonics of the 
permanent current and conductance holds for a much 
larger class of systems. To check this, we have also 
considered Euclidean lattices in section IV Fl namely the 
square and the triangular lattices. 



The application of the strong disorder RG to three 
and four-site clusters is discussed in Appendix [0 In 
agreement with the preceding argument we find a 0o/2 
periodicity for the three-site cluster with a symmetric 
distribution of and a <po periodicity for the four-site 
cluster with ek = 0. The permanent current is expanded 
according to 

JW = ^f]A n sm(2nnpl. (37) 
<Pa ^ V <PoJ 

The numerical results obtained by iterating the RG 
transformations on the fractals with n — 3,4, and on 
the Sierpinski gasket (see Fig.ll5|). are in agreement with 
the general argument in section IV D II The <j>/(j>o har- 
monics vanishes on the fractal with n = 3 and on the 
Sierpinski gasket for a symmetric distribution of on-site 
energies. We obtain a strong 0o harmonics for the fractal 
with n = 4. 

It is interesting to compare these results with stud- 
ies of persistent current in mesoscopic rings. For a sin- 
gle ring this current is expected to oscillate with a pe- 
riod d>n |2flj| . and this prediction has been confirmed in 
two experiments |2ll |22| |. Upon ensemble-averaging with 
a fixed particle number, the periodicity is expected to 
change to tpo /2 j2^J 0, , in agreement with several 
observations on systems with a large number of discon- 
nected rings [2a. 1271 128| | . Ensemble-averaging with a fixed 
chemical potential instead of a fixed particle-number pre- 
serves a dominant 0o periodicity 29] . 

The argument of section lV D ll suggests that this quali- 
tative difference between bipartite and non-bipartite lat- 



E. Two-terminal current 

The evolution of the weight of the 4>q and <fio/2 har- 
monics as a function of e max is shown in Fig. 1161 for the 
fractals with n = 3 and n = 4. We focus here on the 
case r = 0, corresponding to a voltage equal to the max- 
imal coupling. The total current, obtained by evaluating 
the correlation function C' 2 (R) (see Eqs. I|25|) and i|26|l ). 
is expanded in a Fourier series: 

G{4>) = e -rY B n cos 27rnf- . (38) 

The variations of i3i(e max ) and £?2(emax) are shown on 
Fig. 1161 The two-terminal current is almost 4>q /2-periodic 
for the fractal with n = 3 and almost 0o _ P er iodic for the 
fractal with n = 4, which is compatible with the previ- 
ous behavior for the persistent current in section IV PI 
The negative value of B 2 for the fractal with n = 3 im- 
plies a positive magnetoconductance at small flux, which 
looks like weak localization. We deduce from the so- 
lution for a single triangular plaquette in Appendix [0 
that the negative value of B 2 for e max = for the frac- 
tal with n = 3 is due to a i-transformation followed by 
an e-transformation. In this case the </>o-periodicity origi- 
nates from the interference between the two time reversed 
paths around an elementary triangle, reminiscent of weak 
localization. We note also that because of infinite disor- 
der, the largest contribution to the energy comes from 
the first RG transformations that couple to the magnetic 
flux. The permanent current is thus dominated by the 
plaquettes on the smallest scale. 
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n=3 



CO 




FIG. 16: Total current (in units of e 2 /h) evaluated at T — 0, 
averaged over R between R — 15 and R = 20, as a function 
of /, for the fractals with n = 3 and n = 4. The different 
curves within each panel correspond to the ±e model with 
a symmetric distribution of (□ and ■), to the +e model 
with tfc > (o and •), and to the — e model with < (A, 
A). The open symbols in red correspond to the amplitude Bi 
of the 4>q harmonics. The filled symbols in green correspond 
to the amplitude Bi of the </>o/2 harmonics. Note the weak 
localization-like positive magnetoconductance for the fractal 
with n — 3. 



F. Comparison with Euclidean lattices 

Using Green's function methods we calculate the con- 
ductance of an array of nodes interconnected by single 
channel contacts (see Fig. I17|) . Disorder is introduced in 



„(°) 



■Sxk and 



the coordinates (xk,yk) of the nodes: Xk = x k -toxk 

Vk = J/[°' > + $Vk i where the (x^ , ) correspond to the 
coordinates of the nodes of the regular lattice having a 
lattice parameter ao- The bare Green's function depends 
on the distance between the nodes. The variables Sxk and 
Syk are random, and chosen to be uniformly distributed 
in the interval [— W, W], with W < ao/2. A similar model 
has also been investigated in Ref . [33| . Some technical de- 
tails regarding the Green's function treatment are given 
in Appendix iDl The evolution of A\/Aq and A^JAq as 
a function of W/clq are shown on Fig. ^| for the square 
and triangular lattices. The triangular lattice is obtained 
by adding one diagonal on the plaquettes of the square 
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FIG. 17: Schematic representation of the N x N lattice of 
nodes and links used in the Green's function method (with 
N = 4). The lattice is connected to a right and left electrode 
by extended contacts. 



lattice. We obtain a strong </>o/2 harmonics for the trian- 
gular lattice while the two harmonics have approximately 
the same magnitude for the square lattice (see Fig. I18|l . 

For the tight-binding model on the triangular lattice, 
the argument given in section IV Dl can be adapted to 
prove that the Kubo conductivity is </>o/2 periodic, when 
the on-site disorder has an even distribution, and when 
the chemical potential vanishes. In our simulations there 
is an extra phase kpR in the Green's function between 
two nodes (see Ea. lD2|) so that the model that we consider 
is not a tight-binding model. Nevertheless we find a be- 
havior in a qualitative agreement with the general argu- 
ment given in section IVD II the 4>o/2 harmonics is small 
for the triangular lattice, but not vanishingly small. We 
carried out a simulation with g^j — inpo instead of (ID2t , 
corresponding to a half-filled model with kpRk.i = 7r/2, 
and found also a small cf>o harmonics in this case. The 
amplitude of the <f>o/2 harmonics is positive for the tri- 
angular lattice, as opposed to the negative sign obtained 
for the fractal with n — 3 on Fig. This indicates 
that the weak localization-like behavior (associated to 
infinite disorder) breaks down for the Euclidean, conven- 
tional disorder triangular lattice. Similarly, the sign of 
the harmonics Ai of the permanent current is positive 
for the Sierpinski gasket (corresponding to conventional 
disorder) while it is negative for the fractal with n = 3 
(corresponding to infinite disorder). We thus see that the 
magnetoresistive effects in infinite disorder systems are of 
a weak localization-like type while a magnetoresistance 
opposite the one of weak localization can be obtained for 
finite disorder systems. 

The case of the square lattice is qualitatively compati- 
ble with the results obtained by various groups on related 
models. In particular, quasi one-dimensional systems 
have been studied by numerical |30l l3l| or analvtical|32j 
methods. A transition from 0o to (f>o/2 periodicity has 
been found when disorder becomes large enough so that 
the elastic mean free path is comparable to the lattice 
spacing. Similar results for a two-dimensional square lat- 
tice have also been reported 33] . In the strong disorder 
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FIG. 18: Evolution of relative amplitude B\/Bo ■ red) and B2/B0 (•, green) of the 4>o and 4>o/2 harmonics as a function of 
the strength of positional disorder W/ao, for the square (a) and triangular lattices (b). We use a lattice of size N x N, with 
N — 5, and kpao = 5. B2 is positive for the triangular lattice, unlike the case of the fractal with n = 3 on Fig. 1161 



regime, a very interesting 4>o to (b o/ 2 transition has been 
predicted by Nguyen et al. |34tl35| . in a model where local 
random potentials can take only two large and opposite 
values, with probabilities x and 1 — x. At small x, it 
exhibits ^-periodic oscillations in a Aharonov-Bohm ge- 
ometry, and <po /2 appears for a concentration x estimated 
around 0.05 for a square lattice. In all these systems, the 
4>q /2 periodicity occurs when disorder is large enough 
so that propagation paths interfering around a plaquette 
can be assumed to have a random relative phase. 



VI. CONCLUSIONS 

In this paper disordered tight-binding models coupled 
to an external magnetic field were studied on differ- 
ent fractal lattices. The fractals have different type of 
topology (Sierpinski gasket and fractals based on poly- 
gons with n sides), symmetry (bipartite and non bipar- 
tite) and fractal dimension, varying between df = 1 and 
df = 2. We studied the low-energy spectrum of the sys- 
tems, the localization properties (fraction of frozen sites), 
correlation functions and correlation length and trans- 
port properties (permanent current and two-terminal 
current). The main method of calculation is the numer- 
ical implementation of a strong disorder RG approach. 
We generalized the method to non-bipartite lattices and 
to extrinsic diagonal disorder. The results of numerical 
calculations are analyzed by using exact arguments, phe- 
nomenological scaling considerations and by comparing 
with results obtained by Green's function methods on 
Euclidean lattices. To have a further comparison and 
to obtain a classification of random fixed points in frac- 
tals we have also studied the random antifcrromagnetic 
Heisenberg model, as well as the spin-glass model (i.e. 
with random ferro- and antifcrromagnetic couplings) on 
the same lattices and by the same type of strong disorder 



RG method. 

Generally infinite disorder fixed points are observed on 
the fractals based on polygons with n sides, both for the 
random antifcrromagnetic Heisenberg model and for the 
random tight-binding model without extrinsic diagonal 
disorder. On the other hand a strong disorder (i.e. con- 
ventional disorder) fixed point is found on the Sierpinski 
gasket and for the spin-glass (± J) Heisenberg model for 
all types of fractals. The specific properties of low-energy 
fixed points realized in fractals are found in some respect 
different from that already known in ID and 2D regular 
lattices. 

For the random tight-binding model an interesting new 
feature of our study is to consider non-bipartite lattices 
as well as extrinsic diagonal disorder. The strong disor- 
der RG method is shown to generate frozen sites (which 
do not contribute to the average long-range correlations 
and to the conductance) and to create a strongly inter- 
connected effective cluster, which, however is less com- 
pact than that for D > 2 Euclidean lattices. The cor- 
relation function, in the absence of diagonal disorder, is 
shown to exhibit a non-trivial algebraic decay, whereas 
for extrinsic diagonal disorder there is a finite correla- 
tion length, the value of which is related to the cut-off 
of diagonal disorder. The persistent current is shown to 
exhibit a periodicity of (bo/ 2 or </>o, depending on whether 
the elementary plaquettes in the lattices have an odd or 
even number of sites, respectively. The same type of 
periodicity is also observed in the two-terminal current. 
Note that two distinct effects are involved: one is specific 
to bipartite lattices for which only t-transformations are 
carried out. The other is specific to a class of non bi- 
partite lattices for which both e- and ^-transformations 
are carried out but the <pa harmonics vanishes for a sym- 
metric distribution of on-site energies. Lattices with an 
odd number of sites around each plaquette belong to this 
second class of non bipartite lattices, that contains also 
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FIG. 19: (a) RG transformation in which site 1 is frozen, (b) 
RG transformation in which the bond 1-2 is frozen. 

other lattices such as the fractal with n = 3. 

Finally we note that fractals are realized in different 
forms in nature 36] , for instance at a percolation tran- 
sition point. Our results might have some relevance of 
understanding the low-temperature dynamics and trans- 
port in these systems. 



where £1,2 abd £2,1 are complex numbers such that £2.1 = 
(£1.2)* (see Eq. 0}. Considering now a three-site cluster 
{1, 2, 3}, we obtain 

*2,3 = *2,3 — *2,i — £i,3< (A2) 

2. t-transformation 

Let us now suppose that the strongest gap is due to a 
hopping integral £1,2 between sites "1" and "2" and first 
consider the three site clusters {1,2,3}. The renormal- 
ization of £3 is given by 

~ 1 1 

e 3 — e 3 ~ *3.1T *2,3 — *3,2" *1,3' (A3) 

12,1 £l,2 

Considering now a cluster made of sites {1,2,3,4}, the 
renormalization of the hopping amplitude is given by 

1 1 , 

*3.4 — £3,4 — *3.1T £2,4 — £3,2- £1,4- (A4) 

£2,1 £l,2 

3. Derivation of the RG equations 



Acknowledgments 

The authors acknowledge many fruitful discussions 
with J.C. Angles d'Auriac. This work has been supported 
by the French-Hungarian cooperation program Balaton 
(Ministere des Affaires Etrangeres - OM), the Hungar- 
ian National Research Fund under grant No OTKA 
T034183, T037323, T048721, M045596 and M36803. 



APPENDIX A: RG TRANSFORMATIONS OF 
THE TIGHT-BINDING MODEL 

In this Appendix we provide a short derivation of the 
RG transformations of the random tight-binding model 
on a general lattice, not necessarily bipartite. The 
strongest gap can be due either to an on-site energy or 
to a hopping energy. In the former case a fermion or a 
hole is frozen at a given site. In the latter case a fermion 
is frozen in a dimer. 



e-transformation 



As an example we provide a derivation of t- 
transformation in a three-site cluster. The Hamiltonian 
is ho = ti,2C^C2 + t2,i c 2 c i+ e 3 c t c 3' an d the perturbation 
h! is made of the £1,3 and £2,3 terms. The unperturbed 
state is 



1 



1,2 



l*i 



1,0, 1> - |0, 1, 1) 



(A5) 



where the basis is formed by the eight states \n\, ri2, TI3}, 
with 711,^2,713 € {0,1}. The perturbed eigenstate takes 
the form \ip) = \i/jq) + P\ip), where P is the projector in 
the orthogonal to IV'o)- The eigenstate can be expanded 
according to 



\%} + (E-h )- 1 Ph'\iP ) + 



(A6) 



and the correction to the energy is given by AE = 
(t[>o\h' \ip) . After a straightforward calculation we obtain 



AE= - 



1*1,3 



1*2,; 



2ti- 



2*1, 



1 



2U 



-*2,3 — £3,2 



1 



2£i, 



-t 



1,3, 



(A7) 

where the last two terms contribute to the renormaliza- 
tion of e 3 in Eq. (JXSJ). 



Let us suppose that the strongest gap is due to the on- 
site energy ei of site "1" (see Fig. ElK a ))- Considering 
first a two-site cluster {1,2}, the renormalization of the 
on-site energies is given by 

€ 2 = € 2 - £2,1— £1.2, (Al) 
ei 



APPENDIX B: TRANSPORT FORMULA 

In this Appendix we provide a derivation of the trans- 
port formula (12511 . At an energy scale fi, the random 
system consists of "dimers" due to the ^transformations 
coexisting with fermions frozen in a single site due to 
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There is however an additional "local" current: 




FIG. 20: Schematic representation of the system at energy 
scale Q. Dimers of variable length (in red) are formed by the 
t-transformation. The sites not involved in a dimer are either 
frozen by an e-transformation or have not been yet decimated. 
We have specialized two sites a and b in the random system, 
and their corresponding sites a and (5 in the left and right 
electrodes. The sites a and b are connected by a hopping 
t a ,b (dashed line), calculated from the RG. We have repre- 
sented the hopping amplitudes t making an extended contact 
between the random system and the left and right electrodes. 
For clarity we have represented an Euclidian 2D random sys- 
tem but we calculate the conductance in fractal geometries in 
Fig. Q R is the distance between the left and right electrodes. 



the e-transformations. Some of the dimers are frozen 
since they were coupled by \tk,i\ > Other dimers with 
\tk,i | = ^ are being processed by the RG and the remain- 
ing dimers with \tkj\ < f2 correspond to the "active" 
sites 

We consider two 3D normal electrode connected at the 
left and right to the random system by two extended 
contacts (see Fig. 12011 . Two arbitrary sites a and b at the 
two interfaces can be paired by a hopping t a ^ formed 
by a ^transformation. We have t a j, = if a bond a-b 
has not been formed. We note t a ^ a and % a the hoppings 
between sites a and a, and b and (3 (we use t a a = t^ p = t. 
Within perturbation in t a a and tb,p the total differential 
conductance is given by (i) the sum over all pairs of sites a 
and b of the "elastic cotunneling" current I a ,fc associated 
to the hoppings t a ,b, an d (ii) the sum over all sites a of the 
"local" currents T a , a . We expand the elastic cotunneling 
current J a b according to 



1 a ,b - 47T 2 ' rt 2 aja tl f pp a ,oip0,i3 



(Bl) 



where /i?,i(w) = /t(cVr,l — ui) is the Fermi distribution 
function in the right and left electrodes, p a ^ a and pp.p are 
the density of state of the 3D left and right electrodes, 
g^ b and g b R a are the advanced and retarded propagators 
from a to b and b to a, u> is the energy. We suppose that 
voltages Vjr = — V and Vr — V are applied on the left 
and right electrodes. We recognize in Eq. (jB 1|> a term 
proportional to the integral over energy of C' 2 {R) (see 
Eq.EJ. 



3-a,a — 4"7T ^ ^a,aP a ^ a 



+ oo 
-oo 

(B2) 

where fv a ,v b ( LU ) is the out-of-equilibrium distribution 
function in the disordered system, that depends on the 
value of the voltages V a and V&. The distribution function 
fv a y b {oj), determined in such a way as to verify current 
conservation, is in general not equal to the Fermi distri- 
bution. The local current vanishes if 14 = or tb,p = 0, 
because in these cases the out-of-equilibrium distribution 
function fv a ,v b (u) is equal to /l(w). Within the strong 
disorder RG hypothesis the non percolating renormal- 
ized hoppings starting at one interface and not ending at 
the other interface do not contribute to the current since 
the fermions cannot hop from one dimer to the other, 
or from one dimer to one active site. The density of 
states p a ,a in Eq. I|B2|) is thus evaluated only for the per- 
colating bonds frozen at energy f2 while C' 2 (R) involves 
a summation over all percolating bonds that have been 
frozen at energy scales larger than (see sectionlnnnj. 
The "elastic cotunneling" term (|Blll proportional to the 
correlation function C' 2 {R) is thus much larger than the 
"local" term <|ij2t if t = t a . a = tb,p remains finite while 
the area of the contacts is sent to infinity. We have thus 
provided a derivation of the transport formula given by 
Eq. (23. 



APPENDIX C: TIGHT-BINDING MODEL ON 
SMALL CLUSTERS 

Because of strong disorder the dominant contributions 
to the permanent current and conductance are due to the 
RG transformations on the smallest length scale. The pe- 
riodicity of the permanent current can thus be qualita- 
tively understood with simple networks made of a single 
plaquette. 



1. Triangular plaquette 

a. e-transformation followed by a t -transformation 

Considering a triangular plaquette made of three sites 
{1, 2, 3} and pierced by a magnetic flux (/>, we first sup- 
pose that ei is the strongest energy. After eliminating 
the site "1" , the renormalized couplings are given by 

(CI) 

(C2) 

(C3) 

and the correction to the ground state energy after the 
first transformation is S = —lei 1/2. The correction to 



~(1) 

e 2 


= £2 - 


*2,1 


1 

—h,2 

£l 




= £3 ~ 


*3,1 


1 

— *1,3 
£l 




= *2,3 


~t 2 


1 

1 — tl,3 
£l 
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the ground state energy due to the next t-transformation APPENDIX D: GREEN'S FUNCTIONS METHOD 



is 5E^ = -I4 1 !!, with 



A cos (2iT(j)/(po), 

ei 



(C4) 



with A = (4°>) 2 + (t^/exf and B = M^tg, 

where tf 3 ^ corresponds to the value of the hopping in 
the absence of magnetic flux. Assuming a symmetric 
distribution P(ei) = P(— ei), the average ground zero 
temperature state energy given by 



U(4>) = - 



J &iP(ei)[|^( ei )| + $3(-ei)| 



is (/) /2-periodic. 



(C5) 



In this appendix we provide a brief description of the 
Green's function method by which we calculate the con- 
ductance. We introduce nodes at the vertices of a square 
lattice, interconnected by links (see Fig. I17fl . We obtain 
a triangular lattice by adding a link on one diagonal of 
the plaquettes. A node correspond to a Green's func- 
tion gk,k = inpo and a link between the nodes rik and n; 
corresponds to a Green's function 



.'//,'(•»= at] ex p(^/ r ' A - dr ) ' ( D1 ) 



where A is the vector potential and gk,i(4>) is the Green's 
function in the absence of an applied magnetic flux. We 
choose 



b. t -transformation followed by an ^-transformation 

We suppose now that the sites {1,2} are first frozen in 
a dimer. The ground state energy is given by 



^)=-^-I|e 3 -/(^o)| 



(C6) 



with 



t (°) t (°) 

/(0/0 o ) = 2^^ cos (C7) 

Assuming a symmetric distribution of on-site energies, 
we obtain 



U{<t>) 



-t 



(0) 



n,2 ~ l e 3 P(e 3 )de 3 (C8) 

!/«>/« I 
r\fW<h)\ 

\fWM\ / P(e 3 )de 3 . (C9) 

Jo 



The average ground state energy is 0o/2-periodic, and 
the sign of the <j)o/2 harmonics is in agreement with 
Fig.H3(a) fore max = 0. 



2. Square plaquette 

Only i-transformations are performed for a bipartite 
lattice in the absence of extrinsic diagonal disorder. It 
is easy to show that two consecutive ^-transformations 
generate a 0o _ periodic ground state energy. 



9k,i = u p exp ( lk FRk,l), 



(D2) 



where Rk,i is the distance between the nodes rik and ni- 
Eq. (|D2|I corresponds to the Green's function of a bulk 
metal with a Fermi wave- vector kp. We obtained similar 
results with other spatial variations of the bare Green's 
functions. 

The fully dressed Green' s function G is obtained by 
interconnecting the nodes and links by a hopping am- 
plitude T such that nTpo = 1 corresponding to highly 
transparent interfaces. The Green's function G of the 
lattice of connected nodes and links is obtained through 
the Dyson equation G = g + gt&T<®G, where (g> is a short 
notation for a convolution over time variables (becoming 
a simple product after a Fourier transform to the fre- 
quencies), and a summation over the labels of the graph 
of nodes and links. The symbol E denotes the self-energy 
corresponding to the connections between the nodes and 
links. By using two times the Dyson equation we ob- 
tain a linear set of equations by which we determine the 
Green's function Gk,i connecting the nodes and n\. 
The Kubo conductance corresponding to two extended 
tunnel contacts at the left and right sides of the network 
is then proportional to 



5(0) 



(D3) 



k.l 



where k and I run over all sites in the left-most and right- 
most columns respectively. 
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